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Abstract 

We analyze the problem of the Nielsen-Olesen unstable modes in the SU(2) lattice 
gauge theory by means of a recently introduced gauge-invariant effective action. We 
perform numerical simulations in the case of a constant Abelian chromomagnetic 
field. We find that for lattice sizes above a certain critical length the density of 
effective action shows a behaviour compatible with the presence of the unstable 
modes. We put out a possible relation between the dynamics of the unstable modes 
and the confinement. 



1 Introduction 

It is well known since long time that for non-Abelian gauge theories without 
matter fields in the one-loop approximation states with an external constant 
Abelian chromomagnetic field lies below the perturbative ground states. In 
the case of the SU(2) gauge theory the vacuum energy density turns out to 
be [1] in the one-loop approximation 

<H) = + ^,g'H^ In(^) + 0{g^H^) , (1) 

where H is the strength of the external magnetic field. The energy density 
Eq.(l) has a negative minimum at some non zero value of the external field 
H. However, N. K. Nielsen and P. Olesen [2] pointed out that the Savvidy 
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states are not stable due to long-range modes: the Nielsen-Olesen unstable 
modes. 

In a series of papers the proposal was advanced that the stabilization of the 
Nielsen-Olesen unstable modes results in the formation of a quantum liquid, 
the so-called Copenhagen vacuum [3]. 

The problem of the Nielsen-Olesen unstable modes has been reconsidered by 
using variational techniques on a class of gauge-invariant Gaussian wave func- 
tionals [4,5]. It turned out that the unstable modes contribute to the vacuum 
energy density with a negative classical term. In particular, in Ref. [5] it was 
shown that the Nielsen-Olesen modes get stabilized by considering configura- 
tions which differ from the classical external field only in the unstable sector. 
Moreover, the configurations which minimize the energy density screen almost 
completely the classical external chromomagnetic field and contribute to the 
energy density with a negative classical term that cancels the positive clas- 
sical term in Eq.(l). However, it was pointed out that the calculation of the 
one-loop contribution to the vacuum energy density is truly non-perturbative 
and needs a completely non-perturbative approach. 

Several authors studied gauge theories with an external background field on 
the lattice, both in three [6,7] and four [8,9] space-time dimensions. In par- 
ticular, the results of Refs. [6] and [9] showed some evidence indicating the 
presence of the unstable modes. In those papers, however, the external back- 
ground field was introduced via an external source. As a consequence the 
lattice effective action loses the manifest gauge invariance, for the source term 
in the lattice action becomes invariant only under a subgroup of the gauge 
group. 

The aim of the present paper is to analyze the problem of the unstable modes in 
the case of the SU(2) lattice gauge theory by means of the recently introduced 
gauge invariant effective action [10]. 



2 SU(2) in a Constant Background Field 



Let us consider a static external background field A'^^{x). We define the ef- 
fective action which is invariant against gauge transformation of the external 
field by using the so-called Schrodinger functional: 

^{Ar\eM-HT)r\AT), (2) 
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where V projects onto the physical states. On the lattice we can rewrite ^[A^''*] 
as follows 



Z 



A 



cxt 



VU exp(— 5*1^) , 



(3) 



where Sw is the standard Wilson action, and we integrate over the lattice 
links U^{x) with the constraints 



^M(^)U4=o = t/r(^'0) 



(4) 



The external links are related to the continuum gauge field A^^^ by the well- 
known relation 



C/^^t [x) = P exp I Mag £ dt A^^ {x + atfi)^ 



(5) 



with An 



and A^^Hf) 



A'^^{x)Xa/2, and P is the path-ordering 



operator. It should be emphasized that the lattice Schrodinger functional has 
a well-defined continuum limit and, due to the gauge-invariance, does not 
require extra counterterms [11]. The lattice effective action for the external 
background field A^^{x) is given by 
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1 , ZlA"""' 



(6) 



where T is the extension in the Euclidean time. It is straightforward to show 
that, in the hmit T — > oo, r[/l'^''*] reduces to the vacuum energy in presence of 
the given background field. Thus, our effective action is relevant to investigate 
non perturbatively the properties of the quantum vacuum. 

In this paper we are interested in the case of a constant abelian chromomag- 
netic field for the SU(2) gauge theory 



(7) 



Note that, due to the gauge invariance, r[A'^^*] depends on the constant field 
strength Fj^ = H. Thus, the relevant quantity is the density of the effective 
action: 



4 ext 
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In 



(8) 



where Q = V-T, and V is the spatial volume. The external links corresponding 
to A^^ are easily evaluated from Eq.(5). We get 
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^r(f,0) = cos(^)+i.3s.n(^) 

Ur\x, 0) = Ul'^'ix) = Ur\x) = 1 . (9) 

By imposing the periodic boundary conditions we obtain the quantization 
condition 

^ = (10) 
2 Li 

where n^^* is an integer, Li the lattice extension in the Xi direction (in lattice 
units) . 



3 The Unstable Modes 



Let us now briefly discuss the origin of the unstable modes both in the con- 
tinuum and on the lattice. 

In order to evaluate the one-loop effective action in the continuum one writes 

A;{x)^Al{x)+ri;{x) (11) 

where A^{x) — Sfa^^-^XiH and ri{x) is the quantum fluctuation on the back- 
ground field. In the background gauge 

{5-%-ge-^-Al)rf^{x) = Q, (12) 

we rewrite the pure gauge action in the one-loop approximation as 

Sy-m = -S-eiass + ljd'^ r}l{x)Of,ri{x) . (13) 

The one-loop effective action can be obtained by performing the Gaussian 
functional integrations over the quantum fluctuations and including the Faddeev- 
Popov one-loop determinant. However, if we solve the eigenvalue equations 

OtM^) = Ml{x). (14) 
then we flnd that there are negative eigenvalues: 

\u = pl + pl-gH. (15) 

Indeed A„ < when gH > Po + p^- The modes with eigenvalue A„ arc the 
Nielsen-Olesen unstable modes. If we perform formally the Gaussian func- 
tional integration, then the one-loop effective action picks out an imaginary 
part. The point is that in the functional integration over the unstable modes 
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one must include the positive quartic term. If we do that, we find that the 
unstable modes behave like a two-dimensional tachyonic charged scalar field. 
Thus the stabilization of the unstable modes resemble the (dynamical) Higgs 
mechanism. The negative condensation energy we gain in the stabilization 
procedure cancels the positive classical magnetic energy [5]. 

On the lattice we would like to evaluate the Schrodinger functional Eq.(3) 
in the weak coupling region. To this end, we write [11] the lattice version of 
Eq.(ll): 

U,{x) = e^^{iagq^{x))UJ^\x), (16) 
where the fluctuations q^x{x) — q'^{x)a"'/2 satisfy the boundary conditions 

q,{x)U,=o = ^. (17) 

Inserting Eq.(16) into the Wilson action we obtain in the one-loop approxi- 
mation 

Sw = -S^^* + -5^2) (18) 

where 



9^ 



1 — cos 



(f 



(19) 



and S**^^-* is quadratic in the quantum fluctuations q^{x). We are interested 
in the spectrum of the lattice version of the operator in Eq.(14). Un- 
fortunately the discrete eigenvalue equation cannot be solved (analytically) in 
closed form. Obviously it is possible to solve the eigenvalue equation by means 
of numerical methods. We plan to discuss this matter in a future paper. Nev- 
ertheless, if we discard a lot of irrelevant operators in S^'^\ i.e. terms which 
vanish in the naive continuum limit, we flnd an approximate formula for the 
eigenvalues. In particular the unstable mode eigenvalue turns out to be 

A„ = (1 - COSP4) + {1- cos Pa) - sin • (20) 

By using Eq.(lO) and = 27rn^/L^, we get 

A„ = 1 - cos (^^^4^ -M - cos (^^^3^ ~ sin (^^^''''*^ • (21) 



4 Monte Ccirlo Simulations 



We performed our Monte Carlo simulations on lattices with linear sizes Li ~ 
-L4 = 32, L2 — L3 — L and rf^^ — 1. Inserting these values into Eq.(21), it fol- 



5 



lows that Xu < when L > Lait, Lent — 10. Remarkably, these considerations 
suggest that we can switch on and off the unstable modes by varying L. 

We turn now on discussing the outcomes of the numerical simulations. We 
are interested in the vacuum energy density Eq.(8). To avoid the problem 
of computing a partition function, we focus on the derivative of with 
respect to P by taking n''^* (i.e. gH, see Eq.(lO)) fixed. A straightforward 
calculation gives: 



e' 



A 



cxt 



\''^x,n>u'^ /q \"a;,/i>i/^ / ^ext (22) 



where the L/^^(a;)'s are the plaquettes in the (/x, u) plane. As we said, our sim- 
ulations are performed on 32 x x 32 lattices with L = 6, 8, 10, 16, 20, 24. We 
impose periodic boundary conditions. We use standard overrelaxed Metropo- 
lis algorithm to update gauge configurations. The links belonging to the time 
slice ,X4 = are frozen according to Eq.(9). Moreover, we also impose that 
the links at the spatial boundaries are fixed according to Eq.(9). Note that 
this condition corresponds in the continuum to the usual requirement that 
the fluctuations over the background field vanish at the spatial infinity. It is 
obvious that the contributions to ^'[A'^^*] due to the frozen time slice at Xi = 
and to the fixed links at the lattice spatial boundaries must be subtracted, 
i.e. only the dynamical links must be taken into account in evaluating ^'[A*^^*]. 
We denote by fl = L1L2L3L4 the total number of lattice sites (i.e. the lattice 
volume) belonging to the lattice A. Qext are the lattice sites whose hnks are 
fixed according to Eq.(9): 

l^ext = L^L2Ls + (L4 - 1) X (L1L2L3 - (Li - 2)(L2 - 2)(L3 - 2)) . 

(23) 

Hence Qi^t = fl — flext is the volume occupied by the "internal" lattice sites 
(let us denote this ensemble by A). Accordingly we define the derivative of the 
"internal" energy density e{^^ as 



^int 



\i'int ^ / \i'int ^ / fOA\ 



In Figures 1 and 2 we report the derivative of the internal energy density 
versus (3 in units of the derivative of the energy density due to the "external" 
links defined in Eq.(9) 

Several features of Fig. 1 are worth noting. In the strong coupling region 
/3 < 1 the external background field is completely shielded. This behavior is 
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very similar to the one we found for the U(l) gauge theory [10]. Moreover, e[j^^ 
displays a peak at ~ 2.2. These features are displayed also [12] by the specific 
heat in pure SU(2). This is not surprising, for our previous studies [10,13] in 
U(l) showed that e[^^ behaves as the specific heat. 

The derivative of the internal energy density displays the most interesting 
behavior in the weak coupling region j3 > 2.3. Indeed, we feel that this peculiar 
behavior is due to the unstable modes. To see this, we recall that according 
to our previous discussion we have in the weak coupling region 



The first term in Eq.(26) is the classical magnetic energy, while g{rf^^) is the 
one-loop contribution. As a consequence we expect that 



Equation (27) tells us that, in the weak coupling region the derivative of the 

vacuum energy density reduces to Sg^f Indeed, in the U(l) gauge theory this 
is the case [10]. From Figure 1 we see that , for L ~ 6 — 8, e[^^ tends to e^^t 
(in the weak coupling region) as in U(l). However, by increasing L the ratio 
^int/^ext decreases (see Fig. 2). 

In order to ascertain if the effects displayed by our data are not due to the 
lattice artifacts, we contrast in Fig. 3 the derivative of the vacuum energy 
density in the weak coupling region both for SU(2) and U(l) versus the lattice 
linear size L. As we can see, the ratio ^(nt/^ext — ^ it is almost independent 
on L in the case of the C/(l) theory. On the other hand, for the SU(2) case 
Fig. 3 clearly shows that the ratio ^int/^ext — >^ in the infinite volume limit. 
As we shall argue in a moment, this peculiar behaviour can be ascribed to the 
unstable modes. To see this we recall that the stabilization of the unstable 
modes should modify Eq.(26) as follows [5]: 



However, in Ref. [5] it was pointed out that the canceUation of the classical 
energy term happens only in the thermodynamical limit y — > oo. So that we 
expect that Eq.(28) holds on the infinite lattice. On a finite lattice we expect 
that in the weak coupling region 




(26) 




(27) 




(28) 




(29) 
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with X — an/Les, where an — J'^'^/gH is the magnetic length and 



i^eff = ^il: , (30) 

is the lattice effective hnear size. The continuum calculations of Ref. [5] tell us 
that /(O) = 0. In Figure 4 we display the ratio e[^^/e'^^^ for /3 = 5 versus the 
lattice effective linear size. Indeed Fig. 4 shows that e[^^ decreases by increasing 
Z/cff. So that we have e[^^^ when L^s — oo, in agreement with Eq.(28). 
Moreover, we see clearly that a critical length exists such that for L^b < L^q^ 
Eq. (27) holds. It is remarkable that our extimation of the critical length 

LS* « 12 (31) 

is in good agreement with Ref. [9]. Moreover, in the region Les > Lgg* we 
fitted the data trying the power law 

f{x) = kx'^ . (32) 

Indeed we find a rather good fit (x^/d.o.f. ~ 0.3, sohd line in Fig.4) with 

A; = 5(1) , a = 1.43(8). (33) 

In our opinion the existence of a critical length and the cancellation of the 
classical magnetic energy density in the thermodynamical limit constitute a 
rather strong evidence for the unstable modes on the lattice. Remarkably it 
turns out that also the peak value of e{^^ tends towards zero with the same 
power law as implied by Eq. (33). Indeed, in Fig. 5 we show the ratio £(nt/^ext 
taken at the peaks in Figs. 1 and 2. For L^s > -^eff* '^^ ^^'^ (solid line in Fig. 5) 

'-f^ = k'(^y', (34) 

^cxt V-t^cff/ 

where k' = 10(2), a' = 1.5(1) and x^/d.o.f. ~ 0.3. Note that a' = a within our 
statistical uncertainties. Thus Eqs. (29), (32), and (34) imply that £int[^''^*] 
tends uniformly towards zero in the thermodynamical limit. 



5 Conclusions 



Let us conclude by stressing the main results of this paper. 

By using the gauge invariant effective action for background fields on the lat- 
tice we studied non perturbatively the dynamics of the Nielsen-Olesen unstable 
modes. We feel that the existence of a critical length L^g* and the cancellation 
of the classical magnetic energy in the infinite volume limit constitute a strong 
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evidence for the unstable modes on the lattice. Moreover, our numerical re- 



sults are suggesting that e' A^^^ when L^n — >• cxd in the whole /3-region. 
Thus in the continuum limit Leg — oo, /5 ^ cxd the confining vacuum screens 
completely the external chromomagnetic Abelian field. In other words, the 
continuum vacuum behaves as an Abelian magnetic condensate medium in 
accordance with the dual superconductivity scenario. The magnetic conden- 
sate dynamics seems to be closely related to the presence of the Nielsen-Olesen 
unstable modes. In conclusion we feel that our approach to a gauge-invariant 
effective action on the lattice is a useful tool to understand the dynamics of 
color confinement. 



References 



[1] G. K. Savvidy, Phys. Lett. B71 (1977) 133; S. G. Matinyan and G. K. Savvidy, 
Nucl. Phys. B134 (1978) 539; I. A. Balatin, S. G. Matinyan and G. K. Savvidy, 
Sov. J. Nucl. Phys. 26 (1978) 214. 

[2] N.K. Nielsen and P. Olesen, Nucl. Phys. B144 (1978) 376. 

[3] For a review, see: H. B. Nielsen in Particle Physics, eds. I Andric, I. Dadic and 
N. Zovko (North-Holland, Amsterdam, 1981). 

[4] M. Consoli and G. Preparata, Phys. Lett. 154B (1985) 411. 

[5] P. Cea, Phys. Rev. D37 (1988) 1637. 

[6] P. Cea and L. Cosmai, Phys. Rev. D48 (1993) 3364; P. Cea and L. Cosmai, 
[hep-lat 7930600^ ; P. Cea and L. Cosmai, Nucl. Phys. (Proc. Suppl.) 34 (1994) 



234. 

[7] H. D. Trottier and R. M. Woloshyn, Phys. Rev. Lett. 70 (1993) 2053; E. 72 
(1994) 4155. 

[8] J. Ambjorn, V. K. Mitriushkin, and A. M. Zadorozhny, Phys. Lett. B245 (199) 
575; P. Cea and L. Cosmai, Phys. Lett. B264 (1991) 415. 

[9] A. R. Levi, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 161; A. R. Levi and J. 
Polonyi, Phys. Lett. B357 (1995) 186. 



[10] P. Cea, L. Cosmai, and A. D. Polosa, |hep-lat/96010TC 



[11] M. Liischer, R. Narayanan, P. Weisz, and U. Wolff, Nucl. Phys. B384 (1992) 
168; M. Liischer and P. Weisz, Nucl. Phys. B452 (1995) 213. 

[12] B. Lautrup and M. Nauenberg, Phys. Rev. Lett. 45 (1980) 1755; J. Engels and 
T. Scheideler, hep-lat/960741. 



[13] P. Cea, L. Cosmai, and A. D. Polosa, |hep-lat/960702C . 



9 



FIGURE CAPTIONS 

Figure 1. eJjjJA'^''*]/^^^^ versus /9 for L = 6 (stars) and L = 8 (crosses). 
Figure 2. £^int[^'^''*]/^ext versus P for L = 10 (diamonds), L = 16 (triangles) and 
L = 24 (circles). 

Figure 3. The ratio £int[^^^*]/^ext the weak coupling region versus L for U{1) 

(circles, /9 = 3) and SU{2) (squares, = 5). 
Figure 4. The ratio £int[^'''^*]/^ext at /? = 5 versus L^s- The solid line is the fit 

Eq. (32). 

Figure 5. The ratio £^int[^'^''*]/^ext the peak values versus Leg. The solid line is 
the fit Eq. (34). 
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Fig. 1. eJj^^fA^'^^j/egxt versus /3 for L = 6 (stars) and L = 8 (crosses). 



11 



1.75 
S 1-5 

V 

CO 

^^ .25 

0) 

T< 

^ 1 







_ %^ ^ « ^ L=10 



0.5 oo ^ 



A A A A L=16 

O O O O 
0.25 h O L=24 



, A 



XT 



P 



Fig. 2. £jj^^[yl*'^*]/egx^ versus /? for L = 10 (diamonds), L = 16 (triangles) and 
L = 24 (circles). 
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Fig. 3. The ratio e[^^[yl'^^*]/e^j^^ in the weak coupUng region versus L for U (1) (circles, 
/? = 3) and SU{2) (squares, /3 = 5). 
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